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Abstract 

We study the phase transition in generahzed chiral or Stiefel's models using 
Monte Carlo simulations. These models are characterized by a breakdown of 
symmetry 0{N)/0{N — P). We show that the phase transition is clearly first 
order for > 3 when P = N and P = — 1, contrary to predictions based 
on the Renormalization Group in 4 — e expansion but in agreement with a 
recent non perturbative Renormalization Group approach. 

P.A.C.S. numbers: 05.50.+q, 75.10.Hk,05.70.Fh, 64.60.Cn, 75.10.-b 

I. INTRODUCTION 

The critical properties of frustrated spin systems are still under discussion [|l|. In particu- 
lar no consensus exists about the nature of the phase transition in canted magnetic systems. 
One example is the stacked triangular lattice with the nearest neighbor antiferromagnetic 
interactions (STA) with vector spins 0{N) where is the number of spin components which 
is always controversial The non-collinear ground state due to the frustration leads to 

a breakdown of symmetry (BS) from 0{N) in the high temperature to 0{N — 2) in the low 
temperature. This is different from ferromagnets in which the ground state is collinear and 
the BS is 0{N)/0{N — 1). Based on the concept of universality, the class of the transi- 
tion would be different in the two models. We generalize this chiral model for a BS of the 
type 0{N)/0{N — P). We obtain the STA model for P = 2 while we obtain new BS for 
N > P > 3. For example, the case N = P = 3 should correspond to real experimental sys- 
tems, this is also applicable in spin glasses where some disorder is present. Several authors 
have already studied these generalized chiral models applying the Renormalization Group 
technic In mean field ^ for N > P the model shows a usual second order type, but 

for N = P the transition shows a special behavior. This last result could be interpreted with 
the BS in this case being Z2 ® SO{N) and the coupling between the two symmetries leading 
to some special behavior (for example, the case N = P = 2 in two dimensions (d) is always 
very debated |]^). The d = 4 — e expansion gives more information. The picture is very 
similar for all > P > 2 (for details see and references therein). At the lowest order in 
e, there are up to four fixed points, depending on the values of A^ and P. Amongst them are 
the trivial Gaussian fixed point and the standard isotropic 0{NP) Heisenberg fixed point. 
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These two fixed points are unstable. In addition, a pair of new fixed points, one stable and 
the other unstable, appear if the case is > Nc{d) with ^ 



Nc{d) = 5P + 2 + 2a/6(P + 2)(P - 1 



25P2 + 22P - 32 

5P + 2 + 



2J6(P + 2)(P-1) 



(1) 



For P = 2 we find the standard result A^^^ = 21.8 — 23.4 e. On the other hand, for P = 3 we 



obtain = 32.5 - 33. 7e and for P = 4 we obtain A^^ = 42.8 - 43.9 e. A "tricritical" line 
exists which divides a second order region for low d and large from a first order region 
for large d and small A^. From these results Kawamura, using e = 1 {d = 3), obtained that 
Nc{d = 3) < for all P. Thus he concluded that the experimental or numerical accessible 
systems (A^ > 0) were in the second order region. Unfortunately it has been proved that 
the results of 4 — e are, at best, asymptotic |^. They have to be resummed to obtain 
reliable results. Indeed for P = 2 the calculation of the next order in e, combined with a 
resummation technic, leads the experimental accessible systems for A^ = 2 or A^ = 3 in the 
first order region P,p!0|. We believe that the same applies for P > 3. In order to verify 



our assumption, we have done some simulations for P = 3 and P = 4 with N = P and 
A^ = P + 1. The most interesting case is P = A^ = 3 with some possible experimental 
realizations and connection with the spin glasses. Moreover it is meaningful to study the 
generalized model in order to have a better overview. The system we analyze is the Stiefel 



model |]TT|. This model is constructed to have the needed BS. It is closely connected to real 
systems with complicating interactions, which are characterized by the same BS (for the 
case P = 2 see |^ and reference therein). From the principle of universality, models with 
the same BS should belong to the same universality class. Moreover we have shown that 
the use of the Stiefel model allows us to avoid problems which are seen in standard models, 
such as the presence of a complex fixed point (or minimum in the flow) 0,^. 

In the following section II the studied models are presented, we describe the details of 
the simulations and the finite size scaling analysis. Results will be given in section III and 
the last section is devoted to the conclusion. 



II. STIEFEL'S MODELS, MONTE CARLO SIMULATIONS AND FIRST ORDER 

TRANSITIONS 

In this section we introduce different models studied in this work. 

First the V^^^ model which is represented in Fig. |l|. The energy of the model is 

H = JJ:^:[^k{^)■e,{J)] (2) 

ij k=l 

where the P mutual orthogonal A^ component unit vectors efc(z) at lattice site i interact with 
the next P vectors at the neighboring sites j. The interaction constant is here negative to 
favor alignment of the vectors at different sites. Taking a strict orthogonality between the 
vectors is similar to removing " irrelevant" modes corresponding to the variation between the 
spins inside the cell. For example, in the case of a triangular lattice with antiferromagnetic 
interactions (STA) we force the three spins of each cell to have a rigidity constraint with 
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the sum of all the spins being always zero. The obtained model is equivalent to the STA 
at the critical temperature and can easily be transformed into the Stiefel's ¥3^2 model (for 
more detail see 0,0). 

We did not use the clusters algorithm 0] because it gives worse results than the standard 
Metropolis algorithm for first order transitions. 

The method for choosing the random vector depends on the number of components A'^. 
For iV = 3 we follow the method explained in |0] for the direct-trihedral model. We construct 
two orthogonal vectors ei and 62, and the third vector is constructed by the vector product 
of the first two: 

eg = cr ei X 62 (3) 

where a is a random Ising variable, corresponding to the Ising symmetry present in the Vs^a 
model. This is the difference with the direct-trihedral model defined in [0], where no Ising 
symmetry is present. 

To simulate the V^^^ model we follow a similar procedure. We construct now three 
orthogonal unit vectors, = (e^, e^, e^, e^) with k = 1, 2, 3, randomly in four dimensions 
using six Euler angles. The first Oq must be chosen with probability sin(^o)^ d9o, two other 
with probability sin(6'i.2) d6i^2 and the rest three 0^^^^^ with probability ^3,4.5. We obtain for 
ei: 

el = — cos{9^) * cos(6'i) * sin(6'5) * cos(^^2) * cos(6'4) 

— cos(6'3) * cos(6'i) * cos(6'5) * sin(6'4) 

-|- cos(6'3) * sin(6'i) * sin(6'5) * sin(6'2) * cos(6'o) 
-I- sin(6'3) * sin(6'5) * cos(6'2) * sin(6'4) 
— sin(6'3) * cos(6'5) * cos(^^4) 
el = — sin(6'3) * cos(6'i) * sin(6'5) * cos(6'2) * cos(^4) 

— sin(6'3) * cos(6'i) * cos(^^5) * sin(6'4) 

+ sin(6'3) * sin(6'i) * sin(6'5) * sin(6'2) * cos(6'o) 

— cos(6'3) * sin(6'5) * cos(6'2) * sin(6'4) 

-|- cos(6'3) * cos(6'5) * cos(^^4) 
el = — sin(^5) * sin(6'2) * sin(6'o) 

= sin(6'i) * sin(6'5) * cos(6'2) * cos(^4) 

-I- sin(6'i) * cos(^^5) * sin(6'4) 
-I- cos(^^i) * sin(6'5) * sin(6'2) * cos(6'o) 

for 62: 

el = sin(6'o) * sin(^^i) * cos(6'3) 
el = sin(6'o) * sin(6'i) * sin(6'3) 
el = cos(6'o) 
62 = sin(6'o) * cos(6'i) 

and for 63: 
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el = cos(6'3) * sin(6'2) * cos(6'4) * cos(0i) 
+ cos(6'3) * cos(6'2) * cos(6'o) * sin(6'i) 
— sin(6'2) * sin(6'4) * sin(^^3) 
63 = sin(6'3) * sm(6'2) * cos(6'4) * cos(^i) 
+ sin(^^3) * cos(6'2) * cos(6'o) * sin(6'i) 
+ sin(6'2) * sin(^^4) * cos(6'3) 
63 = — cos(6'2) * sin(6'o) 

63 = — sin(6'2) * 003(^^4) * sin(6'i) 
+ cos(6'2) * cos(6'o) * cos(6'i) . 

From this model we can easily create the direct-quadrihedral model and the V4^4, these 
two models being composed of four orthogonal vectors with four components. The differences 
between the two models are the presence of an Ising variables in the V4,4 where right handed 
and left handed are allowed while only one possibility exists in the direct-quadrihedral. The 
direct-quadrihedral and the V4,3 models are topologically equivalent, they should have the 
same low energy physics and therefore belong to the same universality class. The connection 
between the direct-trihedral and the V3^2 models ^ is very similar to the above. We form a 
fourth vector from the vector product of ei, 62 and 63: 

64 = ei X 62 X es (4) 

for the direct-quadrihedral model, and we add a random Ising variable a to the V4^4 model: 

64 = cr ei X 62 X es . (5) 

We follow the standard Metropolis algorithm to update one P-hedral after the other. In 
each simulation between 20 000 to 100 000 Monte Carlo steps are made for equilibration 
and averages. Cubic systems of linear dimensions from L = 10 to L = 25 are simulated. 

The order parameter M for this model is 

^ = ^1:1^^1 (6) 

where Mi is the total magnetization given by the sum of the vectors over all sites and 
is the total number of sites. 

For N = P we define a chirality order parameter: 

^ = j-,e^-{Y[e,) (7) 

where means the vector product x. 

We use the histogram MC technique developed by Ferrenberg and Swendsen which 
is very useful for identifying a first order transition. 

The finite size scaling (FSS) for a first order transition has been extensively studied 
13HT5[. A first order transition can be identified by some properties and in particular by 



the following: 
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a. The histogram P{E) has a double peak. 

b. The magnetization, the chirahty and the energy have hysteresis. 

The double peak in P(E) means that at least two states with different energies coexist in 
the system at one temperature. 

III. RESULTS 

We now present our results for the different models. The Vs^s and the V^^^ show a strong 
first order transition. The hysteresis in E and < M > are shown in Fig. ^ and ^ for the 
Va^a model, and in E and k for the V4^4 model in Fig. ^ and ^ This is in accordance 



with the negative r] exponent for the V^^s model found in [|ri| which describes a first order 



transition because t] must be positive P JT6|JT7[] . This result is understandable because there 
is a coupling between the Ising symmetry and the SO{N) symmetry. We notice that the ¥2^2 
model is also of first order and that the models have a stronger first order transition if N 
is greater (we obtain the same hysteresis for L = 20 for the Vs^s as for L = 10 for the V4,4). 
Thus we can generalize our result that the transition is always of first order for N = P. 

The V4^3 model shows no hysteresis. However a double peak structure appears in the 
energy histogram and becomes more apparent when the size increases (Fig. |^). For greater 
sizes the two peaks are well separated by a region of zero probability, the transition time 
from one state to the other grows exponentially with the size of the lattice. We should 
obtains hysteresis in the thermodynamic quantities when the simulation is not too long. 
The V4^3 model has a first order transition but weaker than the direct-quadrihedral model 
which, for similar sizes, shows hysteresis. As explained above the two models belong to the 
same universality class, i.e. a first order transition, similar to the dihedral model ¥3^2 and 
the direct-trihedral model |0. The addition of the fourth leg to the ¥4^3 model allows the 
first order behavior to be more clearly visible. In Fig. ^ we have plotted our hypothesis for 
the RG diagram flow. Following the initial point, the flow could be under the influence of a 
"complex" fixed point (or minimum of the flow 0) and the system mimics a second order 
transition. Well outside the influence of this fixed point the transition is strongly of first 
order and in the crossover between these two regions the transition is weakly first order. For 
a more developed discussion see 0. 



IV. CONCLUSION 

We have tried to give a general picture of the transition with a 0{N)/0{N — P) break- 
down of symmetry. We have shown by numerical simulations that for N = P = 3 and 
N = P = 4, the transition is clearly of first order. We have generalized our result for all 
N = P. A similar conclusion is obtained for = 4 and P = 3. Using the fact that for 

= 3 and P = 2 the transition is also of first order, we can generalize our result for all 
P = A^ — 1. This is in contradiction with to the conclusion of Kawamura [^] which is based 
on two loops of a 4 — e expansion. As we have noted the e expansion has to be resummed 
to obtain reliable results. We can try to achieve this by forming simple Fade approximants. 
For a function / = a + 6e we obtain the approximation / = 1/(1 — efe/a) which we apply 
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to e = 1 (rf = 3) and P =2, 3 and 4. We obtain Nc{P = 2) ~ 10, N^{P = 3) ~ 16 and 
Nc{P = 2) ~ 21. Unfortunately the results can not be perfect and in particular the result 
for P = 2 is not close enough to the result including the next order expansion N^. = 3.39 
10| which demonstrates that the lower-order e expansions are useless in this case. However 



we remark that the Nc "resummed" increases with P which is in agreement with our result, 
i.e. that the initial point in the renormalization flow is farther away from the mimic of the 
second order region Thus the systems will show a stronger first order transition. This 
result matches with a recent study on the case P = N = 3 which is based on a non pertur- 
bative Renormalization Group procedure |TB|. We conclude that transitions for N = P and 
N = P + 1 are of first order for all A^. 
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FIGURES 
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FIG. 1. Stiefel's model Vs^s and their inter- 
actions. ei(i) interacts only with ei(j), neither 
with 62 (j) which interacts with 62 (i), nor with 
63 (i) which interacts with 63(2). 



FIG. 3. Magnetization versus T for the 1/3^3 
model. The size of the system is L = 20. See 
comments in Fig. |2| 




FIG. 2. Internal energy per spin E versus 
T for the V3^3 model. Lines are guides to the 
eye. The arrows indicate if the MC simulation is 
cooling (circle) or heating (square) the system. 
The system size is L = 20. A hysteresis is clearly 
visible. 



FIG. 4. Internal energy per spin E versus 
T for the 14,4 model. The system size is L = 10. 
See comments in Fig. ^ 
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FIG. 5. Chirality versus T for the ¥4^4 FIG. 7. Internal energy per spin E versus T 
model. The system size is L = 10. See com- for the direct-quadrihedral model. The system 
ments in Fig. size is L = 20. See comments in Fig. |2| 
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FIG. 6. Energy histogram P{E) as a func- 
tion of the energy per site E for the V^^^ model 
for various sizes L at different temperatures of 
simulation Tl: Tig = 1.1802, T20 = 1.1771, 

= 1.1758. 



FIG. 8. Hypotheses of Hamiltonian flows 
induced by renormalization-group transforma- 
tions. The arrows (dashed line) indicate the di- 
rection of flow under iteration. 
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